Magnet field sensing beyond the standard quantum limit under the effect of decoherence 
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Entangled states can potentially be used to outperform the standard quantum limit which every classical 
sensor is bounded by. However, entangled states are very susceptible to decoherence, and so it is not clear 
whether one can really create a superior sensor to classical technology via a quantum strategy which is subject 
to the effect of realistic noise. This paper presents an investigation of how a quantum sensor composed of many 
spins is affected by independent dephasing. We adopt general noise models including non-Markovian effects, 
and in these noise models the performance of the sensor depends crucially on the exposure time of the sensor 
to the field. We have found that, by choosing an appropriate exposure time within non-Markovian time region, 
an entangled sensor does actually beat the standard quantum limit. Since independent dephasing is one of the 
most typical sources of noise in many systems, our results suggest a practical and scalable approach to beating 
the standard quantum limit. 



Entanglement has proven itself to be one of the most in- 
triguing aspects of quantum mechanics, and its study has lead 
to profound advances in our understanding of physics. Aside 
from these conceptual advances, the exploitation of entangle- 
ment has lead to a number of technical advances both in com- 
putation and communication 1 1 1, and more recently in metrol- 
ogy |2|. In quantum metrology, entanglement has been used 
to demonstrate enhanced accuracy both in detecting the phase 
induced by unknown optical elements and for accurately mea- 
suring an unknown magnetic field. It is this latter case which 
is the focus of the present paper, and so we will adopt the 
terminology of field sensing. 

In order to estimate an unknown field, one usually prepares 
a probe system composed of L distinct local subsystems, ex- 
poses this to the field for a certain time, and measures the 
probe. Comparing the input with the output of the probe gives 
us an estimate of the field. Importantly, there is an uncertainty 
in the estimation, and this uncertainty is related to how the 
probe system is prepared. When the probe system is prepared 
in a separable state, the uncertainty decreases as l/\/Z EHU 
by the central limit theorem, a scaling known as the standard 
quantum limit. On the other hand, by preparing a highly en- 
tangled state, it is in principal possible to achieve an uncer- 
tainty that scales as 1 /L, known the Heisenberg limit |I2]|5J. 

Most of the literature on quantum sensing focuses on using 
photons to probe an unknown optical element. Using a NOON 
state ||6]-(8], it is possible to measure the unknown phase shift 
with higher resolution than the standard quantum limit. An 
L-photon NOON state can achieve a phase L times as large as 
than in the case of a single photon over the same channel. Re- 
cent publications f9 - 12] have considered an analogous tech- 
nique for field sensing with spins. In this paper, we consider 
an experiment involving a probe consisting of spin-| systems. 
The spin qubits can couple to the magnetic field and therefore 
one can estimate the value of the magnetic field by the probe. 
In order to obtain higher resolution than the standard quantum 
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limit, one can use a GHZ state, as has been demonstrated in 
recent experiments llOUTTl [T3l . 

In solid state systems, one of the main barriers to realising 
such sensors is decoherence, which degrades the quantum co- 
herence of the entangled states. GHZ states, in particular, are 
very susceptible to decoherence, and decohere more rapidly as 
the size of the state increases lfT4l[T5l . Therefore, it is not clear 
whether a quantum strategy can really outperform an optimal 
classical strategy under the effects of a realistic noise source. 
The effect of unknown but static field variations over the L 
spins has been studied by Jones et al ifTOl . The uncertainty 
of the estimation depends on the exposure time of entangled 
states to the field, where they are affected by noise, and they 
have found that, for an optimal exposure time in their model, 
the scaling of the estimated value is i which beats the stan- 
dard quantum limit. However, the underlying assumption that 
the fields are static could be unrealistic for many systems, as 
actual noise in the laboratory may fluctuate with time. Huelga 
et al have included such temporal fluctuations of the field in 
their noise model 1 16] and have shown that GHZ states cannot 
beat the standard quantum limit under the effect of indepen- 
dent dephasing by adopting a Lindblad type master equation 
lfT6l . Even for the optimal exposure time, it was shown that 
the measurement uncertainty of a quantum strategy has the 
same scaling behavior as the standard quantum limit in their 
noise model. Since independent dephasing is the dominant er- 
ror sources in many systems, these results seem to show that, 
practically, it would be impossible to beat the standard quan- 
tum limit with a quantum strategy. 

However, the model adopted by Huelga et al is a Markovian 
master equation |TT| which is valid in limited circumstances. 
The Markovian assumption will be violated when the correla- 
tion time of the noise is longer than the characteristic time of 
the system. For example, although a Markovian master equa- 
tion predicts an exponential decay behavior, it is known that 
unstable systems show a quadratic decay in the time region 
shorter than a correlation time of the noise fT8'/T9l. In this 
paper, we adopt independent dephasing models which include 
non-Markovian effects and we investigate how the uncertainty 
of the estimation is affected by such noise. We have found 
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that, if the exposure time of the entangled state is within the 
non-Markovian region, a quantum strategy can indeed provide 
a scaling advantage over the optimal classical strategy. 

Let us summarize a quantum strategy to obtain the Heisen- 
berg limit in an ideal situation without decoherence. A state 
prepared in |+) = ^(|0) + |1)) will have a phase factor in 
its non-diagonal term through being exposed in a magnetic 
field and so we have ^(|0) + e~'*''|l)) where 5 denotes the 
detuning between the magnetic field and the atomic transi- 
tion. On the other hand, when one prepares a GHZ state 
|V'>GHZ = 73(1'^'^ • • • 0) + |11 • • • 1)) and exposes this state 
to the field for a time t, the phase factor is amplified linearly 
as the size of the state increases as 
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|^(<)) = ^(|00...0)+e-^"|ll...l)) 



(1) 



Therefore, the probability of finding the initial GHZ state after 
a time t is given by P = | + | cos{LtS). In practice, one may 
use control-not operations to map the accumulated phase to a 
single spin for a convenient measurement I.10J . The variance 
of the estimated value is then given by 

P{1~P)/N 



A'5 = 



(2) 



\dP/d6\^ 

where N is the number of experiments performed in this set- 
ting lfT6l . For a given fixed time T, one can perform this ex- 
periment T/t times where t is the exposure time, and so we 
have N = T/t. Hence we obtain |A^(5| — and so the 
uncertainty in 6 scales as L~^, the Heisenberg limit. 

First let us consider the decoherence of a single qubit. 
Later, we will generalize to GHZ states. Our noise model 
represents random classical fields to induce dephasing. We 
consider an interaction Hamiltonian to denote a coupling with 
an environment such as 

Hi = A/(i)<T. (3) 

where /(<) is classical normalized Gaussian noise, (j^ is a 
Pauli operator of the system, and A denotes a coupling con- 
stant. Also, we assume symmetric noise to satisfy f{t) = 
where this over-line denotes the average over the ensemble 
of the noise. When we solve the Schrodinger equation in an 
interaction picture, we obtain the following standard form 

Pi{t)- Pa^^{-i>^)" dti dt2--- I dU 



[Hi{t,),[Hj{t-2 



(4) 



where po = is an initial state and pi{t) is a state in 

the interaction picture. Throughout this paper, we restrict our- 
selves to the case where the system Hamiltonian commutes 
with the operator of noise, as this constitutes purely dephas- 
ing noise. By taking the average over the ensemble of the 
noise, we obtain 

P/W-Po = V^HA)" / dt, dt2--- dt^ 
,^1"! Jo Jo Jo 



where [ctz , [ct^ . • • • , [ct^, po] • ■ ■ ]]n denotes the n-folded com- 
mutator of Po with (7z- Since all higher order cumu- 
lants than the second order are zero for Gaussian noise, 
/(^i)/(^2) • ■ • f{tn) can be represented by a product of cor- 
relation functions 120^. Therefore, the decoherence caused by 
Gaussian noise is characterized by a correlation function of 
the noise. For Markovian white noise, a correlation function 
becomes a delta function while a correlation function becomes 
a constant for non-Markovian noise with an infinite correla- 
tion time such as 1// noise. To include both noise models as 
special cases, we assume the correlation function 
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where Tc denotes the correlation time of the noise. In the limit 
Tc — > 0, this correlation function becomes a delta function, 
while in the limit of Tc ^ 00 it becomes constant. 

Since f {ti) f {12) ■ ■ ■ }{tn) can be represented by a product 
of correlation functions, we obtain 



PI 



<^z, [<Jzi ■■■ , V^z,Po\ ■ ■ ■ \\2n 



n=0 



s.s' = ±l 



where \s) is an eigenvector of (t^. Also, ^{t) denotes the sin- 
gle qubit decoherence rate defined as 



l{t) = 



'(-1 + e -g) 



4A- 



'Tcerf( — ) 



(8) 



where erf(a:) is the error function. Note that, for t ^ Tc, the 
decoherence rate becomes constant as ^{t) ~ -iX^Tc- So, in 
this regime, we can derive a Markovian master equation from 
(|7]i, which has the same form as adopted by Huelga et al fW\. 



dpi{t) 
dt 



-AV, 



(9) 



(5) 



Note that, although our model can be approximated by Marko- 
vian noise in the long time limit (t ^ Tc), we are interested 
in the time periods t ^ Tc and t ^ Tc where non-Markovian 
effects become relevant. 

For an initial state |+) = ^(|0) + |1)), the non- 
diagonal terms of the density matrix show a decay behavior 
of (0|p7(t)|l) = ^e\p[—j{t)t]. For t ^ Tc, the state shows 
an exponential decay (0|p/(t)|l) ~ On the other 

hand, we have (0|p7(t)|l) ~ ^ (quadratic decay be- 
havior) for t <C Tc , which is the typical decay behavior of 
1// noise | 2T|j23l . The behavior of the decoherence rate is 
illustrated in Fig.[T| 

We next consider decoherence of a GHZ state induced by 
random classical fields. Extending the Hamiltonian in ([3]), the 
interaction Hamiltonian denoting random classical fields for a 
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FIG. 1: We have plotted the behavior of a single qubit decoherence 
rate for 1/f noise -^X^t (the green line), Markovian noise 4A Tc 
(the red line), and our classical noise model (blue dotted line) 7(t) 
defined in l[8j where A and Tc denote a coupling constant with the en- 
vironment and a correlation time of the noise, respectively. For an L- 
qubit GHZ state, the decoherence rate becomes L times greater than 
the value plotted here. Note that our noise model shows a transition 
from a non-Markovian quadratic decay to an exponential Markovian 
decay. Here, we fixed parameters as A = 0.25 and Tc = 1. 



many-qubit system is as follows: 

L 



(10) 



1=1 



where fi{t) denotes the noise acting at site I and has the same 
characteristics as /(t) mentioned above. Also, since we con- 
sider independent noise, we have fi{t)fii{t') (x Sij>. When 
we let a GHZ state be exposed to a magnetic field under the 
effect of the random magnetic fields, the state will remain in 
thesubspace spanned by {^^j^ |0);and(^^j^ because we 
are only considering phase noise. Thus, we can use the same 
analysis as for a single qubit. In the the Schrodinger picture, 
we will obtain 



Pit) 



1 



^ »Lt5-L7(t)t((g)|0)j(l|) 



1 = 1 



1=1 



^-iLt5-L'Y{t)t I 



1=1 



(11) 
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where S is the phase induced by the magnetic field to be mea- 
sured. It is worth mentioning that the decoherence rate for 
a i-qubit GHZ state becomes L times the decoherence rate 
7(t) of a single qubit. Therefore, the variance of the estimated 
value becomes 



A^<5 = 



„27(t)Lt 



1 



1 



TLHsin^LtS) TLH 



(12) 



where we use (j2|i, and therefore we obtain the following in- 
equality by using x > sin x > -x for < x < |: 



„27(t)Lt „ 1 1 



TLH 



< 



(13) 



These inequalities depend on both the system size L and 
the choice of exposure time t. We wish to see if there is any 



choice of t for which we beat the standard quantum limit. We 
find that this limit can indeed be beaten provided that we chose 
shorter t values for larger systems. For example, suppose that 
we chose t according to the rule t — sL~^, where s is a con- 
stant with the dimension of time, and z is a non-negative real 
number whose optimal value we will determine. Then, from 
dSll, we have 



l7(t) 



< 



(14) 



for large L and hence the decoherence rate scales as 7(t) = 
Q{L~^). Throughout this paper, for a function f{L), we say 
f{L) = 8(i") if there exists positive constants J and K such 
that JL" < f{L) < KL" is satisfied for all L. 

In ( pjj i, the term in the exponential ^{t)tL goes to infinity 
as L cx) for z < 1/2, and so the uncertainty |A^(5| di- 
verges, which means that a large GHZ state becomes useless 
to estimate 6. Therefore, we consider the case of z > 1/2. By 
performing a Taylor expansion of e^''^*^^*, we obtain 



„27(t)Lt _ I 



E 



{2-f{t)Lty 



n=l 

for z > 1/2 where we use ( p^ . 
So, from ( [T3| l and ( [T5| ), we obtain 



(15) 



(16) 



for z > |. Therefore, when z 



i, we achieve a scaling of 
the uncertainty |A(5| = 9(L^i) and this actually beats the 
standard quantum limit. Note that this is the same scaling as 
the magnetic sensor under the effect of unknown static fields 
studied in lITOl . and so our result for the fluctuating noise with 
time becomes a natural generalization of their work. 

The decoherence model described above is a classical one. 
Next, we make use of a quantized model where the environ- 
ment is modeled as a continuum of field modes. The Hamil- 
tonian of the system and the environment are defined as 



H = 



2 

1=1 



l,k 



E^ 

Lk 



9kCrz{l)h]J)Lk (17) 



and b\ denote annihilation and creation operators 



where hi^k 

for the bosonic field at a site I 



Also, since we consider in- 
dependent noise, we assume that 6; fc commutes with h[, ^, for 
I ^ I'. This model has been solved analytically by Palma et al 
ll24ll and the time evolution of a GHZ state is given by 

= 2((g)|0)/(0|) + 2e'^*'""^<*^*((g)|0>,(l|) 
1=1 1=1 

+ Je^"*'^^r('^*((g)|lMO|) + i((g)|l>,(l|), (18) 



1=1 



1=1 



where T{t) denotes a decoherence rate defined as r(t) = 
J log(l + oj'^t^) + J log C''"^^^^*-* ) where we have taken the 
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Boltzmann constant fc^ = 1. Here, T and ujc denote the 
temperature and cut off frequency respectively. As we take 
/cs = 1, the temperature T has the same dimension as the fre- 
quency LUc- We have a constant decoherence rate T{t) ~ 2tiT 
for t ^ T^^, which signals Markovian exponential decay, 
while we have V{t) — uj^t for t <C w^T^, which is the charac- 
teristic decay of 1// noise. So, by taking these limits, this 
model can also encompass both Markovian noise and 1// 
noise. Given the calculations in the previous section, the vari- 
ance of the estimated value for the magnetic field becomes 



A^5 = 



TLH-sin^{Lt&) ^ tUH 



(19) 



By performing a calculation exactly analogous to the case of 
the random classical field considered above, one can show 
that the scaling law for the uncertainty becomes A^(5 = 
for z > ^ when we take an exposure time t — 
s/L^ . On the other hand, the uncertainty will diverge as L 
increases for z < i. Therefore, by taking z — the uncer- 
tainty scales as 9(i^i) and so one can again beat the stan- 
dard quantum limit as before. 

We now provide an intuitive reason why the uncertainty of 
the estimation diverges for large L when z is below ^ in both 
of noise models. It has been shown that an unstable state al- 
ways shows a quadratic decay behavior in an initial time re- 
gion 1 18 , 19 1, and therefore the scaling behavior of the fidelity 
of a single qubit should be F = = 1-^2+0(^3) 

where C is a constant. So the scaling behavior of the fidelity 
becomes F = 1 — CLt^ + 0{t^) for multipartite entangled 
states under the effect of independent noise |25]. If we take 
a time as < = sL^^, to first order we obtain an infidelity of 
1 — F ~ Cs^L^^^^. So this infidelity becomes larger as L 
increases for z < ^, which means coherence of this state will 
be almost completely destroyed for a large GHZ state. On the 
other hand, as long as we have z > ^, the infidelity can be 
bounded by a constant even for a large L and so the coher- 
ence of the state will be preserved, which can be utilized for a 
quantum magnetic sensor 

Finally, we remark on the prospects for experimental real- 
ization of our model. To experimentally realise such a sensor, 
one has to generate a GHZ state, expose the state in a mag- 
netic field, and measure the state, before the state shows an 
exponential decay. Although it has been shown that an unsta- 
ble system shows a quadratic decay behavior in the initial time 
period shorter than a correlation time of the noise lll8liT9l . it is 
difficult to observe such quadratic decay behavior experimen- 
tally, because the correlation time of the noise is usually much 
shorter than the typical time resolution of a measurement ap- 
paratus for current technology. After showing the quadratic 
decay behavior, unstable systems usually show an exponential 
decay 1 18 1 and, in this exponential decay region, it is not pos- 
sible to beat the standard quantum limit |J_6J. However, it is 
known that 1 // noise has an infinite correlation time and one 
doesn't observe an exponential decay of a system affected by 
1// noise |23 |. Therefore, a system dominated by such noise 



would be suitable for the first experimental demonstration of 
our model. For example, it is known that nuclear spins of 
donor atoms in doped silicon devices, which have been pro- 
posed as qubits for quantum computation |26|, are dephased 
mainly by 1/f noise |!26lj27 | and so they may prove suitable to 
demonstrate our prediction. 

In conclusion, we have shown that, under the effect of inde- 
pendent dephasing, one can obtain a magnetic sensor whose 
uncertainty scales as Q{L^i) and therefore beats the stan- 
dard quantum limit of L^s. We determine that, to outper- 
form a classical strategy, the exposure time of the entangled 
states to the field should be within the non-Markovian time 
region where the decoherence behavior doesn't show expo- 
nential decay. Since the noise models adopted here are quite 
general, our results suggest a scalable method to beat the stan- 
dard quantum limit in a realistic setting. 
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